Terahertz Kerr and Reflectivity Measurements on the Topological Insulator Bi 2 Se 3 
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We report the first terahertz Kerr measurements on bulk crystals of the topological insulator 
Bi2Se3. At T=10K and fields up to 8T, the real and imaginary Kerr angle and reflectance mea- 
surements utilizing both linearly and circularly polarized incident radiation were measured at a 
frequency of 5.24 meV. A single fluid free carrier bulk response can not describe the line-shape. 
Surface states with a small mass and surprisingly large associated spectral weight quantitatively fit 
all data. 
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I. INTRODUCTION 

The discovery^ - — and realization of 3D strong topologi- 
cal insulators (TIs)2r— have generated much recent inter- 
est. A strong spin-orbit coupling can produce a topolog- 
ically protected metallic surface state which has a Dirac- 
like dispersion with a non-degenerate chiral spin state 
that can not be destroyed or gapped by (non-magnetic) 
impurity scattering. Since the spin direction of the quasi- 
particle is locked to the momentum direction^ suppres- 
sion of Anderson localization results from the destruc- 
tive interference associated with the quaisparticle's Berry 
phase upon backscatteringJ^ With the action of local 
ordered magnetic moments, the surface state becomes 
gapped.— In zero magnetic field with the Fermi energy 
within the gap, a quantized Kerr angle associated with 
the half quantized Hall conductance is a predicted signa- 
ture of the topological origin of the surface stated 

These and other predicted properties have stimulated a 
recent flurry of activity to probe the exotic surface state. 
Although the surface states have been observed by angu- 
lar resolved photo-emission spectroscopy (ARPES),— ~— 
unambiguous separation of bulk from surface state con- 
duction has plagued transport characterization of the 
surface state due to the rather large intrinsic charge den- 
sity in the bulk conduction band in as-grown Bi2Te3 and 
Bi2Se3. Various schemes to circumvent this difficulty 
in bulk crystals have either implemented counter dop- 
ing during growth^ minimized Se vacancies^ or fine 
tuned annealing processes^ One intermediate goal has 
been to sufficiently deplete the bulk to allow gating or 
chemical doping of the surface in suitably thin crystals 
to move the chemical potential of the bulk into the gap^ 

There are several advantages of terahertz Kerr mea- 
surements over other transport measurements which al- 
low separation of the surface from the bulk response. The 
surface states are expected to have a a different mass than 
the bulk, depending upon the location of the chemical 
potential, giving rise to two different resonant cyclotron 
frequencies. This effect can be further exaggerated by 
application of an electric field either from gating or chem- 
ically doping the surface causing 'band bending' effects, 
pushing the Dirac point closer to the chemical potential. 



Secondly, the line-shape as a function of magnetic field 
associated with a surface state is qualitatively different 
than the bulk response. Thirdly, no electrical contacts 
are required. Lastly, Kerr measurements on thick crys- 
tals allow probing of a single surface state. 

In this article, Kerr measurements are reported on 
two bulk crystals of Bi2Se3 at an optical frequency of 
5.24 meV in fields up to 8T. Line-fits of the data show 
evidence of surface states. For a measured surface state 
cyclotron mass of m cs = 0.033 m e , the associated spectral 
weight is about a factor of 20 times larger than expected 
based upon the corresponding number density estimated 
from ARPES data. Measurements of other samples from 
the same as-grown crystal show a higher carrier con- 
centration. Alternatively, a simple model of the optical 
response which assumes the presence of two crystallites 
with different carrier concentrations can qualitatively de- 
scribe all of our data. Future work to optically charac- 
terize carrier density homogeneity as well as an analysis 
of surface state signals expected by gating are discussed. 



II. EXPERIMENTAL DESCRIPTION 

A. Kerr System Description and Calibration 

For calibration purposes, the Faraday rotation and 
ellipticity induced by a sample in magnetic field was 
measured utilizing a polarization modulation technique 
with a heterodyne detection scheme, a technique de- 
tailed clsewhereJi The complex Faraday angle and rela- 
tive transmission were concurrently measured. 

The instrument was modified allowing simultaneous 
measurement of the complex Kerr angle and relative re- 
flection by the insertion of a beam splitter, polarizer, and 
collection optics which include an off-axis parabolic alu- 
minum mirror and a 2 K silicon composite bolometer as 
depicted in Fig.Q] The Kerr and Faraday angle are mea- 
sured separately. 

Assuming a rotationally invariant sample, the complex 
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FIG. 1: (Color online) Schematic of the Faraday and Kerr op- 
tical system capable of operating on various laser lines in the 
terahertz spectral region. The figure, adapted from Ref.H7l 
which details the Faraday experiment, depicts additional op- 
tics within the red-dashed trapezoid allowing Kerr measure- 
ments by capturing reflected light from the sample. 



Kerr and Faraday angle are defined as follows: 

tan(6>A') = i(r+ — rJ)/{r+ +r_) 
tan(0 F ) =i(t+-t-)/(t+ + t-) 

where r (t) is the reflection (transmission) complex Fres- 
nel amplitude of the sample expressed in a circular po- 
larization basis. The + and - represent right and left 
circular polarization states, respectively. 

Absolute measurement of the real part of the Faraday 
and Kerr angles is verified by inserting a linear polarizer 
into the system immediately downstream from the rota- 
tor. The measured rotation angle equals the polarizer 
angle for small angles which is a very good approxima- 
tion for the signals measured in Bi2Se3. Large angle cal- 
ibration may be accomplished in a similar manner, but 
corrections from higher order terms involving the Fresnel 
amplitudes must be explicitly taken into account in the 
modulated signal^ 

Calibration of the imaginary part of the Faraday and 
Kerr angles is accomplished by measuring the known re- 
sponse of a GaAs 2-DEG heterostructure. A NiCr antirc- 
flection (AR) coating was deposited on the backside to 
minimize multiple reflections within the GaAs substrate. 
There is a well defined relationship between the real and 
imaginary parts of the Faraday angle as well as the Kerr 
angle. Therefore, the imaginary part may be calibrated 
against the real part. In this way, calibration of the sys- 
tem is performed as a function of magnetic field in the 



same manner as samples of interest without the addition 
of other optical components. Utilizing the formalism de- 
veloped in Sec. lA2l the real and imaginary part of the 
Faraday and Kerr angle data shown in Fig.[5Ja-d) are 
simultaneously fit where the mobility fx, cyclotron mass 
m, 2D number density n, and sheet resistance of the AR 
coating i?rj are free parameters. The B-field value of 
the resonance determines the mass m, the width of the 
resonance determines the scattering rate characterized by 
the mobility fj,, the overall strength of the response deter- 
mines the number density n, and deviations of the AR 
coating sheet resistance Rq away from the ideal value 
tends to mix the real and imaginary parts of the angles 
giving rise to anisotropics in the resonant peaks. 

The fit parameters are found to be the following: 
n=2.0 x 10 11 cm -2 , fi = —1.19 x 10 5 cm 2 /Vsec where 
the negative sign denotes the sign of the carriers, 
m=0.0695m e , and R=1.01x Kar where m e is the mass of 
an electron and Kar is the sheet resistance for the perfect 
AR coating condition of GaAs given by Z /{nc a As — 1) = 
145 fi/D where Zq = 377 f2/D is the impedance of free 
space and ncaAs = 3.6 is the index of refraction for GaAs. 
The mass is consistent with the known mass of GaAs^ 
and the antireficction coating is very close to ideal. The 
fits are overlayed on top of the reflection and transmis- 
sion data shown in Fig.[2je,f). Essentially, all 6 curves 
are very well reproduced with only two parameters, n 
and /j,. Together with the polarizer calibration, the 2- 
DEG GaAs calibration indicates that both the Kerr and 
Faraday experiments perform exceptionally well. 

This particular data set also illustrates the robust- 
ness of the Kerr system. Instability associated with the 
two-laser system occurred between +8 and ~+5.5T ap- 
pearing as a loss of power and commensurate increase 
of relative detector noise in the Kerr and reflection data 
sets. Although the reflection data was visibly affected, 
the Kerr signal was essentially unaffected. The reflec- 
tion measurement is derived by the ratio of the signals 
from two detectors ('Kerr detector' and 'reference detec- 
tor' in Fig.[T]) and therefore susceptible to large spacial 
mode fluctuations, but the Kerr signal is derived from 
the modulated signal measured by a single detector. 



B. Bi2Se3 sample description 

Bulk crystals of Bi2Se3 were grown by varying the heat- 
ing conditions of the melt and stoichiometric ratio of bis- 
muth and selenium*^ Surfaces were cleaved perpendicu- 
lar to he c-axis, and the backsides were cut and ground 
yielding a total bulk thickness of about 0.5 cm. Sam- 
ples were mounted such that the applied magnetic field 
was parallel to the c-axis and the radiation propagation 
direction. The effective area measured was defined by 
the optical spot size ~ jzfl.Omm. During a time span of 
about 10 seconds, the surfaces of the crystals were cleaved 
using scotch tape in air and immediately inserted into a 
vacuum of ~ 10~ 7 Torr. 
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FIG. 2: (color online) Calibration of the Faraday and Kerr angle measurement system. The Faraday angle, Kerr angle, and 
relative transmission and reflection from a GaAs 2-DEG heterostructure were measured at an optical frequency of 5.24 meV 
and temperature of 80 K as a function of magnetic field. The Faraday and Kerr angles are simultaneously fit with effectively 
two parameters, the number density and mobility, found to be n=2.0 x 10 11 cm -2 and /i = —1.19 x 10 5 cm 2 /V sec. The fits (red 
line) are overlayed on top of the data (black dots) for comparison. The same two parameters are used to calculate reflection 
and transmission (normalized to the B=0 value) and overlayed on the data. 



III. KERR AND REFLECTION 

MEASUREMENTS ON Bi 2 Se 3 AND MODEL 
DESCRIPTION 



model with an associated 2-D surface plasma frequency 
Lo ps , surface cyclotron mass m cs , and surface scattering 
rate j s . Details of the model are discussed in Sec. lA 31 



The Kerr angle and concurrent reflection measure- 
ments of linearly polarized light are shown in Fig.|5fa- 
c) and (e-g) for two samples both measured at 10 K in 
fields up to 8 T. A separate measurement was performed 
on Sample 1 where the waveplate in the rotator depicted 
in Fig. [1] was statically set producing elliptically (nearly 
circularly) polarized light. The reflection was then mea- 
sured and normalized to the B=0 T value. The data is 
reported in Fig.[3|d). ft characterizes the degree of mix- 
ing between left and right circular polarization states via 
(1 — /3)R+ + /3R_ where R± is the reflection produced by 
purely circularly polarized light. In this measurement, (5 
is found to be 0.13. 

To analyze the results, a model is constructed where 
the bulk is considered to be a thick slab whose optical 
response is characterized by a dielectric function consist- 
ing of a Drude carrier response with an associated un- 
screened bulk plasma frequency io v b , bulk cyclotron mass 
m c b, and bulk scattering rate 7& along with the addition 
of two TO phonon modes whose parameters have been 
previously determined^ The surface state resides at the 
interface whose admittance is modelled within a Drude 



IV. DISCUSSION 
A. Analysis of the optical response 

Sample 1 data shown in Fig.|3fa-d) can not be quan- 
titatively nor qualitatively reproduced by a homogenous 
bulk response with any arbitrary carrier concentration, 
cyclotron mass, or scattering rate. The blue dash-dot 
curves in Fig.[3][a-d) show a representative example of 
line-shapes expected from a simple bulk response, calcu- 
lated with the specific parameters oo p b = 420 cm -1 (cor- 
responding to a number density of n& ~ 3.2 x 10 17 cm -3 ), 
76 = 12 cm , and m c i, = .135 m e . Although the high 
field Kerr data in Fig.[3ta,b) is well described by this par- 
ticular set of parameters, strong qualitative deviations 
exist at fields below \B\ < 4T suggesting that another 
carrier response with a much smaller cyclotron mass may 
be superimposed. The reflection data in Fig.|3{c,d) also 
show very different behavior at low fields \B\ < 4T com- 
pared to a bulk response. The small humps which ap- 
pear at fields around ~ ±1 T in Fig.^c) and ~ —IT in 
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FIG. 3: (color online) Kerr angle and reflection of linearly polarized light normalized to the B=0T value were measured at 
10 K and an optical frequency of 5.24 meV. Reflection of nearly circularly polarized light normalized to the B=0T value was 
a separate measurement performed only on Sample 1. (a-d) Sample 1 data (black dots) is depicted with several different 
fits superimposed. A homogeneous bulk response can not describe the data (blue dash-dot). A homogeneous bulk response 
with surface states (red solid) fit the data. The weighted optical response of two crystallites with different carrier density 
(green dash-dot) can reproduce the salient features of the data. Slightly elliptically polarized light, characterized by parameter 
/3=0.13, is used in calculations for Fig.|3jd). (e-h) Sample 2 data (black dots) show a very different response but was only 
2 mm away from Sample 1 in the as-grown crystal. A homogeneous bulk response (blue dash-dot) does not quite describe the 
data. Addition of a surface state response does not improve the fits. A weighted response of two crystallites with different 
carrier density (green dash-dot) more closely resembles the data. No reflectivity data using circular polarized light was taken 
for Sample 2. /3 is assumed zero in calculations in Fig.[3{h). 



Fig-EJd) can not be qualitatively reproduced with any 
value of bulk carrier concentration, cyclotron mass, or 
scattering rate. 

The addition of a surface state with parameters uj ps = 
41 cm -1 , 7 S = 28 cm -1 , and m cs = 0.033 m e results in 
the red solid curves in Fig. [3J a-d) which fit all the data 
remarkably well. The surface state cyclotron mass cor- 
responds to a resonance at B=1.5T (for a frequency of 
5.24 meV). The scattering rate associated with the sur- 
face is substantially larger than that of the bulk, a re- 
sult which might be expected based upon Shubnikov dc 
Haas measurements on similarly grown samples^ and 
linewidths estimated from ARPES datajii Reasonable 
linefits may be achieved with a surface state plasma fre- 
quency as low as uj ps = 20 cm -1 although deviations from 
the data begin to become more pronounced. 

From the ARPES measured dispersion of the surface 
states (see Eq. lA4p ^ the cyclotorn mass given by Hkf/vf 
and number density given by 7rfc?/(27r) 2 may be parame- 
terized in terms of the Fermi energy. A mass of 0.033 m e 
corresponds to a Fermi energy of ^20meV above the 
Dirac point, a number density of 10 11 cm -2 , and cor- 
responding surface plasma frequency of 1.3 cm" 1 . The 



bulk plasma frequency corresponds to a Fermi energy of 
approximately 20 meV above the bulk conduction band 
edge which is itself located ~ 200 meV above the Dirac 
point. 

ARPES observations correspond to a large occupation 
number of the surface states, n s ~ 10 12 cm _2 jii This 
amount of charge acting on the bulk carriers would de- 
plete the free carriers to a depth d = n s /ni, ~ 100 nm 
(where n& is the bulk number density). Since ARPES 
concurrently observed the bulk bands but the measure- 
ment technique is sensitive to depths of only ~lnm, it 
must be assumed that there is a fixed neutralizing posi- 
tive surface charge coming from excess Se + or adsorbed 
ions, or possibly an intrinsic effect. Therefore, band 
bending effects are apparently small in vacuum cleaved 
crystals. However, for samples exposed to ambient condi- 
tions, band bending may be more severe^ possibly due to 
oxidation of the surface. This chemical doping moves the 
Dirac point toward the chemical potential and leads to a 
depletion region in the bulk, presumably with d < n s /nb- 
Since the optical penetration depth is much larger than 
100 nm at terahertz frequencies^ an insulating depletion 
layer has no significant optical response and the experi- 



ment sees the surface carriers only as a modification to 
the bulk. 

More significantly, the surface plasma frequency ob- 
tained from the fits to the optical response are about a 
factor of 15 to 30 larger than predicted from conventional 
transport theory based on the ARPES measured surface 
state dispersion relation. A possible explanation is that 
the crystal has multiple surface layers within the opti- 
cal penetration depth, perhaps separated by interstitial 
selenium, dividing the crystal into ~10 thin flakes which 
are optically very close but electronically decoupled. The 
cleave process may favor the region with layers of flakes 
where these interstitial selenium may be more prevalent 
offering a weaker link between planes. 

Given the unconventional spin-chiral nature of the 
eigenstates, which imply spin currents perpendicular to 
the charge currents, an unconventional matrix element 
for Landau level transitions may also lead to enhance- 
ment of the spectral weight, a possibility that needs fur- 
ther investigation. 

As a possible alternative explanation, the structure ob- 
served at low field in the optical response of Sample 1 
may arise from spatial variation of the bulk charge car- 
rier concentration. Consider the possibility that the op- 
tical spot covers a region where two crystallites of dif- 
ferent carrier concentrations exist. Modelling this sce- 
nario, a simple weighted average of two individual bulk 
responses without surface states is performed with the 
following parameters: one bulk response weighted by 0.84 
with parameters uipbi = 334 cm -1 , 7 p {,i = 6.5 cm -1 , and 
fncbi = 0.126 with a second bulk response weighted by (1- 
0.84) with parameters uj p b2 = 548 cm -1 , 762 = 8.2 cm -1 , 
«V>2 = 0.148 m e . The green dash-dot curves in Fig.[3ja- 
d) show all of the salient features of the data which was 
attributed to the surface states. The difference between 
the two bulk masses is unexpected for a simple parabolic 
band near a gapped band edge. The coarse qualitative 
features of the data may be reproduced by choosing the 
bulk masses to be the same with the other parameters 
chosen similarly although the fits arc substantially worse. 
Therefore, even though the fit is not very good quantita- 
tively, it is possible that a distribution of carrier concen- 
trations may be the culprit giving rise to the observed 
line-shapes. 

Perhaps lending more credence to the spatial varia- 
tion of carrier concentration argument is that a com- 
pletely different response was measured by a separate 
sample which was only 2 mm away from Sample 1 in the 
~ 01.5 cm as-grown crystal. The measurement results of 
Sample 2 are shown in Fig.[«Hc-g). A single bulk response 
qualitatively describes all the features of the data, as de- 
picted by the blue dash-dot curves, whose parameters are 
ujpb = 583cm , 76 = 7 cm" 1 , and m C 6 = 0.145 m e . j3 is 
assumed zero in calculations in Fig.[3Jh). This is a much 
higher bulk plasma frequency than found for Sample 1. 
Addition of a surface state response does not improve 
the fits. However, modelling the bulk carrier concentra- 
tion variations as was done for Sample 1 where one bulk 



response is weighted by 0.70 with parameters u> p bi = 
583 cm -1 , 7(,i = 5.3 cm , and m C 6i = 0.145 m e and a 
second bulk response weighted by (1-0.70) with param- 
eters u> p b2 = 250cm , 762 = 15cm , m C 62 = 0.145m e 
yields a good quantitative fit to all three data sets as 
shown by the green dash-dot curves in Fig.[3je-g). A 
third sample measured 2 mm away from Sample 2 (and 
4 mm away from Sample 1) yielded the same response as 
sample 2. 

Controlling carrier concentrations to well within a fac- 
tor of ~ 10 over cm-sized crystals is very challenging. 
Using the Kerr and reflectivity measurements, it is not 
possible to determine whether such a response originates 
from a distribution of carrier concentration in a single 
crystallite, or the presence of two (or more) crystallites 
with different concentrations. 



B. Other methods to characterize surface states 

Using the Kerr measurement technique as a routine di- 
agnostic of sample carrier concentration variation is not 
feasible since it is resource intensive. Since inhomogene- 
ity is an important general issue, a routine diagnostic 
tool is required. One method presently being explored 
utilizes FTIR spectroscopic reflection measurements. As 
representatively shown in the reflectivity spectral mea- 
surements in figure 4 of Ref.QjJ a large dip occurs at 
the plasma edge. The frequency location of the dip gives 
information regarding the carrier concentration. Mea- 
surements on different parts of a sample utilizing a small 
spot size ~ .00.5 mm in the infrared spectral region is cur- 
rently being investigated as a diagnostic tool to quantify 
the degree of carrier concentration variation. 

Another method that may be favorable for observing 
surface states presently being developed utilizes a gate. 
If the surface states have a cyclotron mass and optical 
spectral weight commensurate with the ARPES mea- 
sured dispersion, the line-fitting procedure as used for 
the above Kerr and reflection data for discerning the sur- 
face state will most likely be extremely difficult with the 
present signal-to- noise of the system, even for fairly small 
surface scattering rates (see Sec. IB II for details). 

By gating, the surface state optical response can be 
modulated with negligible effects produced by the bulk. 
Any modulated signal coming from the bulk would nec- 
essarily be due to the screening layer after the depletion 
layer which is of thickness ~ l/k s , the Thomas- Fermi 
screening length. k s ~ 15 nm at n& = 10 17 cm~ 3 , a 
typical carrier concentration. The charge in this layer 
n-d ~ nb/k s ~ 10 11 cm" 2 which is much less than the 
surface state carrier density. More important, however, 
is that k s does not depend on gate voltage and there- 
fore is not modulated. Secondly, by modulating a gate, 
signal-to-noise is expected to increase by over an order of 
magnitude. Averaging between gate-on and gate-off con- 
figurations in fixed field on relatively short time scales 
eliminates magnetic field drift issues and minimizes time 



drift issues. Thirdly, due to this large gain in sensitiv- 
ity, even very high surface scattering rates should be de- 
tectable. Expected signals produced by gating, modelled 
in Sec. lB 21 are expected to be well within experimental 
capability. Fourthly, the expected Kerr response from the 
surface states for bulk plasma frequencies ranging from 
200 to 600 cm -1 have qualitatively the same line-shapes. 
Even with possible carrier inhomogencitics, the detection 
of surface states is possible by modulating a gate with 
modest shifts of the surface state Fermi energy. How- 
ever, characterization of the surface states will require 
reasonably homogeneous carrier concentrations. Lastly, 
developing the technique is a necessary stepping stone for 
future experiments: measurement of the half-quantized 
Hall conductivity in topological insulators inside the gap 
created by a magnetic layer iL 



V. CONCLUSION 

Kerr measurements using heterodyne detection in the 
terahertz regime are reported. Measurements of a 2-DEG 
GaAs heterostructure prove the system works exception- 
ally well. Results on Bi2Se3 crystals show that the tech- 
nique is sensitive to bulk carrier concentration. For one 
particular sample, the Kerr and reflectance data are very 
well described by including surface states. The measured 
surface spectral weight uj ps (with an associated surface 
cyclotron mass m cs = 0.033 m e ) is ^20 times larger than 
expected for a single surface state based upon ARPES 
measurements. Spatial variation of carrier concentration 
can also qualitatively account for the data. Separation 
of bulk and surface states may be accomplished in the 
Kerr geometry with greater sensitivity by incorporating 
a gate, a proposed next step to clearly separate the bulk 
and surface state effects. 
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Using Maxwell's equations for radiation at normal inci- 
dence, boundary conditions on the electric and magnetic 
field may be applied across any boundary between two 
dielectric slabs i and i + 1 of thicknesses dbi and dbi+i 
with bulk complex indices n.; and nj+i and a metallic 
thin film at the interface with a sheet conductivity of 
<7j dfi where dfi is the thickness of the film. The vec- 
tor notation (E r ,E z ) represents right and left travelling 
waves. The electric field may be expressed in either a lin- 
ear polarization basis notated as (E£,E^.) and (ET yl Ej y ) or 
in a circular polarization basis as (E+,E+) and (E!L,E|_)- 
To derive the equations generally, we use the notation 
(E£,E-,) to represent one of the polarization basis states 
y, x, +, or -. 

The matrix Mj iQ which transfers an arbitrary E-field 
(E^,E^) across an interface consistent with Maxwell's 
equations is given by the following: 



M 



~ 2"i+ 



-X 
Tli+l,a ~ 

t- J1i+1,Q 



Ui,a 
VVi.a 






l,a " 



Vi,ot 



where the admittance of the thin film is y^ a = Zo<Ji t0t dfi 
and Zq is the impedance of free space. 

The matrix K ia transfers an arbitrary oscillating E- 
field (Eq,E„) at frequency ui across a dielectric of thick- 
ness dbi and complex index rii jQ is given by the following: 



K 



t l ni a dbi uj/c 







-i^ iQ dbi uj/c 



In general, an incident right travelling wave will be 
assumed incident on the first interface with a resulting 
net reflected travelling wave to the left from the first 
interface and a net right travelling transmitted wave from 
the last interface. The Fresnel transmission and reflection 
amplitudes for a stack of dielectrics of m interfaces with 
thin metallic films are then found by the relation: 



M n 



| Ki+\, a Mi >a 



i=i 



GaAs 2-DEG 



Appendix A: Multilayer response analysis 

1. General Notation 

The formalism for deriving the Fresnel reflection and 
transmission amplitudes for a system consisting of mul- 
tiple layers of 'dielectric' slabs with thin films at the in- 
terfaces are derived for a general number of layers. The 
model will be applied to two specific cases: a GaAs 2- 
DEG heterostructure with a metallic antircflcction coat- 
ing and a bulk Bi2Sc3 crystal with a surface state. 



The above formalism can be applied to calculate the 
Kerr and Faraday angle as well as the reflection and 
transmission for the calibration sample. A GaAs 2-DEG 
is modelled as a thin metallic sheet whose sheet conduc- 
tance is described within a simple Drudc model located 
at the interface between vacuum and a GaAs dielectric 
slab. The backside of the sample consists of an interface 
between a GaAs dielectric slab and vacuum with an an- 
tireflection coating, a thin metallic NiCr film, modelled 
within a Drude model. The antircflcction coating is fully 
characterized by the dc sheet resistance, Rq = 1/(0 df), 
since the scattering rate of NiCr is very large compared 
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FIG. 4: (color online) Calculated optical response produced by the bulk properties and the expected surface state response as 
measured by ARPES of Bi2Se3 as a function of magnetic field for series of surface state Fermi energies measured above the 
Dirac point as given in the legend: E F =0, 10, 20, 30, 40, 50, 75, 100, 150, 200, 250, 300, 350, and 400 meV. The Kerr angle as 
well as the reflectance (normalized to zero field) in a pure circular or linear state are calculated. Each row of graphs represents 
a different set of Drude parameters describing the bulk and surface state response. The second and third row graphs, (e-h) and 
(i-1), are offset for clarity such that the Kerr angle is zero and R/R[0T] is unity in zero magnetic field. The graphs with the 
Fermi energy Ef = correspond to the Dirac point where there is no contribution from the surface state, only bulk response. 
All graphs are calculated with a bulk mass m c (, = 0.15 m e and a bulk scattering rate 7b = 10 cm -1 and using the bulk phonon 
parameters as measured in Ref.fig. (a-d) Very low bulk density to demonstrate the response expected from the surface states: 
io p b = 10cm -1 , 7 S = 10cm -1 . (e-h) to p b = 400cm~ 1 , 7 3 = 3cm _1 (i-1) u p b — 400 cm -1 , 7 S = lOcm -1 . (m-p) to p b = 400cm~ 1 , 
7 S = 25 cm -1 . 



to the optical frequency. The circular polarization basis 
is easiest to work with in applied magnetic field. 

The parameters are as follows: ni,± = 1, n.2,± = 3.6, 
7i 3 ,± = 1, w = 5.24 meV, (Z adf)i,± = Wp al /(7i - i(u> ± 
u! cs i)), db 2 = 0.5mm, and {Zq<j df)2,± = Zo/R\j where 
Wcsi = eB/{mc) is the cyclotron frequency of the 2-DEG, 
7i = e/(m/Lt) is the scattering rate of the 2-DEG, uj ps i = 
47rne 2 /(mc) is the plasma frequency of the 2-DEG, m is 
the quasiparticle mass, B is the applied magnetic field, c 
is the electronic charge, c is the speed of light, /j, is the 



mobility, and n is the 2-D number density. 

The Fresnel amplitudes r± and t± as a function of n, /i, 
Rrj, m, and B are then found by substituting the above 
parameters into the following equations: 



Y) = M 2 , + K 2>+ M h+ (7 
f ) = M 2 ,-K 2 ,-Mi,. ' 



(Al) 



The reflection R and transmission T for unpolarizcd 



light, reflection R± and transmission T± for circularly 
polarized light, and the Faraday 6p and Kerr 6x angles 
are calculated from the Fresnel amplitudes as follows: 



R ± = r±r*± 
T± = t±t* ± 
R = (1/2) (R + - 
T = (1/2) (T+ 4 
tan(9p) = i{t+ - 
tan(#x) = i(r+ 



■R_) 

T-) 

■t-)/(t+ + t-) 
-r-)/(r+ + r-) 



(A2) 



where '*' notates the complex conjugate. 



3. Bi 2 Se 3 

Bi2Sc3 is modelled as a surface state Drude response 
on top of a very thick bulk slab whose response consists 
of a free carrier Drude term with two phonon resonances. 
The bulk phonon parameters as well as the approximate 
cyclotron mass and scattering rate have previously been 
measured.— 

The bulk dielectric function is given by the following: 



U)(U1 ± W cb - I'Jb) 



E 



Sj 



uj 2 - uj(uj - I'jj) 



(A3) 
where j enumerates the phonon oscillators, cjj is the TO 
phonon frequency, Sj is the phonon spectral weight, a 
subscript b refers to bulk electrons, e^, is the dielec- 
tric constant at higher infrared frequencies, uj p b is the 
unscreened plasma frequency for free electrons: u^ = 
inN^e 2 /m^, where iVj, e and m^ are concentration, 
charge, and an effective mass of electrons; 7^ and jj 
are the damping rates; u c b = eB/(m c bc) is the bulk cy- 
clotron frequency, where B is magnetic field, m c b is bulk 
cyclotron mass, and c is the speed of light. 

The reported phonon parameters are uj\ = 67 cm -1 , 
Si = 3.94 x 10 5 cm" 2 , 71 = 5cm" 1 , cu 2 = 134cm" 1 , 
5 2 = 7056 cm" 2 , 72 = 2 cm" 1 , and eoo = 25.6.i2 The 
reported bulk cyclotron mass is approximately m c b w 
0.16 m e where m e is the bare electron mass and the ap- 
proximate bulk scattering rate 7^ ~ 8 cm" 1 , although 
both are left as parameters in the model. The relation 
ml = m c b is assumed true, so that Nj, (or equivalently 
the plasma frequency, uj p b) is the third adjustable pa- 
rameter associated with the bulk response. With these 
values, e± is a function of radiation frequency u> and ap- 
plied magnetic field B, and the bulk plasma frequency 
ujpb, mass m c b, and scattering rate 7;,. 

Since the surface state penetration depth is ~nm 
which is much less than the optical wavelength for ter- 
ahertz radiation, the surface state may be modelled as 
a two dimensional sheet of charge. The model of a 
1 cm thick Bi2Sc3 slab in vacuum with a Drude sur- 
face state response at the first interface are as follows: 
7H,± = 1, n 2 ,± = a/c±, «3,± = 1, (^0O"4f)l,± = 
u ls/(ls -i(u±u cs )), db 2 = 1cm, (Z adf) 2 ,± = where 



uj cs = eB/(m s c) and u ps = (4:ir/c)n s e 2 /m s . Even for 
very low densities of bulk carriers, the properties of the 
last interface do not effect the reflected optical response 
due to the large attenuation within the thick bulk crystal. 

Inserting the above values into Eqs. lAll give the reflec- 
tion and transmission Fresnel amplitudes as a function 
of w, B, uj p b, 76, m cb , uj ps , 7 S , and m cs . The Fresnel 
amplitudes arc also given in Rcf.|l9|, Eqs. 4 to 6. 

Digitizing and fitting the reported ARPES measured 
dispersion relation for the surface state from Ref.Hll gives 
the following functional form: 



E F = 2.02k f + 10.44/c? 



/ 



(A4) 



where the Fermi energy Ep is measured in eV above the 
Dirac point and kf is measured in A~ l . n s (or w ps ), m cs , 
Vf , and kf are fully determined by the dispersion relation 
for any given Fermi energy. 

Using the dispersion relation, R, T, R±, T±, 8p, and 
9k are calculated from Eq. lA2l as a function of u, B, w p b, 
76, m c6 , Ef, and j s . 



Appendix B: Expected optical signatures of surface 
states based upon ARPES measurements 

1. Expected Kerr and Reflectance signals 

Figure [4] shows calculated results for the expected Kerr 
and reflection signals, both from linearly and circularly 
polarized light, as a function of magnetic field for vari- 
ous Fermi energies ranging from to 400 meV above the 
Dirac point. All graphs in Fig. [4] are calculated assuming 
a bulk cyclotron mass value of m c b = 0.15 m e and bulk 
scattering rate jb = 10 cm" 1 . The surface state plasma 
frequency uj ps and cyclotron mass m cs are determined by 
the ARPES measured dispersion (Eq. lA4[) for a specific 
surface Fermi energy. Each row of graphs are calculated 
with a different bulk plasma frequency u p b and surface 
scattering rate 7 S . 

The optical response for a pure surface state with 
7 S = 10 cm" 1 is shown in Fig.QJa-d) (the response of 
the bulk is negligible for cu p b = 10 cm" 1 ). All signals 
are well within the capability of the Kerr measurements. 
Note that Rc(6k) arid Im(#Rr) show antiresonant- and 
resonant-like behavior. The cyclotron mass and strength 
of the response increase with Fermi energy. The reflection 
in purely circularly polarized light depicted in Fig.[4jd) 
exhibits a hump in the cyclotron resonant active mode. 
The reflection in linear polarization depicted in Fig.|^c) 
(composed of equal amounts of left and right circularly 
polarized light) shows two humps since there is an ac- 
tive cyclotron mode in positive and negative fields. The 
qualitative structure of the Kerr and reflection signals 
is instructive since they appear superimposed on larger 
bulk responses at larger values of ujpb- 

Row 2 of figures, namely Fig.Q|e-h), show the optical 
response expected with a bulk plasma frequency of uj p b = 
400 cm" 1 and surface scattering rate 7 S = 3 cm" 1 . The 
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FIG. 5: (color online) Calculated optical response produced by a hypothetical gate voltage which shifts the Fermi energy of 
the surface state by 30meV. The change in Kerr angle produced by this shift at various Fermi energies for a specific surface 
scattering rate is shown for each column. All graphs are calculated using the bulk phonon parameters measured in Ref.QiJ as 
well as the Drude bulk parameters to p t — 400 cm" 1 , yt, — 10 cm" 1 , and m c j, = 0.15 m e . The surface state 2-D plasma frequency 
and mass are calculated based upon the Fermi level and the measured ARPES dispersion relation.— The surface scattering 
rates are assumed to be (a)7 s = 10cm -1 , (b) 7 S = 25cm -1 , (c) 7 S = 50cm -1 , and (d) 7 S = 100cm -1 . 



bulk is depicted by the black curves since no carriers are 
in the surface state at the Dirac point, Ef = 0. The 
curves are offset for clarity: in zero magnetic field, the 
Kerr angle is always zero and the normalized reflection is 
always equal to unity. With such a low surface scattering 
rate, the sharpness of the surface state response is seen 
clearly superimposed upon the bulk response. 

The next row of figures, Fig-HJi-1) , show the effect of 
increasing the surface scattering rate to 7 S = 10cm -1 . 
Again, the curves are offset for clarity. Although the ef- 
fect of surface states is discernable in many of the curves, 
fitting the lineshapes with the current signal-to-noise of 
the measurement system becomes problematic. 

The last row of figures, Fig.@|m-p), show the effect of 
increasing j s to 25 cm -1 . Fitting lineshapes on this level 
to discern the effects produced by surface states is not 
conducive to the current measurement scheme. 



2. Expected Kerr signals with application of a gate 

A gate which concurrently allows ample throughput 
optical power while inducing negligible Kerr or Faraday 
rotation has been accomplished. High resistivity silicon 
substrates with a 300 nm film of Si02 were ion implanted 
with peak concentrations > 10 19 /cm 3 resulting in a den- 
sity profile thickness of ~ 0.3 /im. Characterization of 
these devices show very small Faraday rotation up to 8 T, 



excellent optical throughput power, no carrier freeze out, 
and large break down voltages. The devices are currently 
being used to gate large CVD grown graphene in similar 
magneto-optical experiments. 

Estimates of the shift in Fermi energy of the surface 
state associated with the application of electric fields 
commensurate with the typical breakdown voltages of 
Si02 depends upon the Fermi energy. For currently 
available samples the surface state charge density is suffi- 
ciently large that it dominates the charge transfer to the 
sample due to the gate. As a result the band offset is rea- 
sonably given by the change in the Fermi energy due to 
the charge transfer to the surface states. The depiction 
depth in the bulk is then given to a good approxima- 
tion by d 2 = 2Uoeo/e 2 rib (where Uq is the energy above 
the Dirac point and eo is the bulk static dielectric con- 
stant). Therefore, the net charge which one can add or 
subtract from the surface state without breakdown of the 
insulator is very nearly equal to the maximum allowable 
charge transfer with the application of a gate. Using 
a maximum applicable electric field value of 3.3 x 10 8 
V/m for SiC>2 is equivalent to a net change in surface 
charge density of 7 x 10 12 cm -2 . Based upon the ARPES 
measured dispersion, this gives a shift of Fermi energy 
AE f «90meV at E f =0, AE f «30meV at E f =200meV, 
and AEf «15meV at .E/=350mcV. Since the gate volt- 
age can be reverse biased as well, roughly twice these 
changes can be realized. The resulting depiction depth 



10 



decouples the surface states, but is small enough that the 
insulating depletion layer is optically invisible at THz fre- 
quencies. 

Based upon these estimates, the application of a gate 
is chosen to produce a change in Fermi energy of A.Ef = 
30meV. Fig. [5] shows the expected Kerr signals predicted 
by the ARPES measured surface state properties with 
bulk parameters of u p b = 400 cm -1 , jf, = 10 cm -1 , and 
m c & = 0.15 m e . Each graph shows Fermi level shifts 
of 0-30meV, 30-60meV, 70-100meV, 120-150meV, 170- 



200 meV, 270-300 meV, and 370-400 meV. Graphs (a-d) 
show the effects of increasing the surface scattering rate 
from 10 cm" 1 to 100 cm -1 . Reflectance measurements 
are not very sensitive, showing that R/R[0T] <0.5% for 
7s = 25 cm" 1 and <0.1% for 7s = 50 cm- 1 . 

Changes in Kerr angle at the level of O.lmrad is very 
reasonably obtainable in this configuration. It should 
be noted that the signals are maintained at these gen- 
eral levels even for higher bulk plasma frequencies up to 
600 cm- 1 . 
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